Visualizing the Spin of Individual Molecules 
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Low-temperature spin-polarized scanning tunneling microscopy is employed to study spin trans- 
port across single Cobalt-Phathalocyanine molecules adsorbed on well characterized magnetic 
nanoleads. A spin-polarized electronic resonance is identified over the center of the molecule and ex- 
ploited to spatially resolve stationary spin states. These states reflect two molecular spin orientations 
and, as established by density functional calculations, originate from a ferromagnetic molecule-lead 
superexchange interaction mediated by the organic ligands. 



PACS numbers: 72.25.-b,73.20.At,75.70.Rf,85.65.+h 

Conceptually new device structures accounting for siz- 
able quantum effects will be needed if the downscaling of 
electronic and magnetic devices were to continue. One 
of these new concepts is the marriage of molecular elec- 
tronics and spintronics, where functional molecules be- 
come active device components within a circuitry where 
information is carried by spins [TJ [2]. Progress toward 
this tantalizing goal rely on our understanding of spin 
transport and magnetism in reduced dimensions, where 
fundamental playgrounds extend to the extreme limit of 
single atoms and molecules. Studies include spin trans- 
port across a well chosen molecule sandwiched between 
two magnetic leads [3j [4], or even atomic-size constric- 
tions formed by bringing two leads into contact [5j EJ. 
One of the limiting drawbacks is the variability of the 
resulting conductance, which comes from the incomplete 
knowledge we have of the molecule- lead interface. For in- 
stance, only a few experimental studies have focused on 
the interaction between a molecular spin and a magnetic 
substrate [7J [8]. A better understanding would enable 
us to target the chemical engineering needed for building 
the desired spintronic functionalities into a molecule. 
In the past years, it became possible with spin-polarized 
(SP) scanning tunneling microscopy and spectroscopy 
(STM and STS) to directly observe the interplay be- 
tween magnetism and surface structure with atomic res- 
olution. SP-STM can also serve as a model tunneling 
magnetoresistance device since the junction includes two 
well defined magnetic leads - the tip and the sample - 
separated by a vacuum barrier. A link can then be accu- 
rately established between spin transport and density of 
states. Recently, through SP-STM it was possible to ev- 
idence how the magnetization switching of a nanocluster 
is influenced by the spatial location of the spin injection 
[9]. Another example can be found in SP-STM of single 
atoms [10]. When the SP current flows across a mag- 
netic atom adsorbed on a magnetic surface rather than 
directly into the surface, the tunneling spin transport is 
significantly affected, and some control can be exerted 
through the choice of the atom. 

Here we show how tunneling spin transport can be 



modified by "dressing" atoms with organic ligands. We 
combine low-temperature SP-STM and model calcula- 
tions to study a model system consisting of individual 
Cobalt-Phthalocyanine (CoPc, Fig. [T^i) molecules ad- 
sorbed on a magnetic substrate. The interaction between 
the molecule and the substrate results in two molecu- 
lar spin orientations clearly discernible in the differential 
conductance of the SP tunnel junction (dl/dV). Our 
results provide a direct visualization of stationary spin 
states, that differs from the molecular Kondo state where 
spin flips occur in time [TTJ [I2j [I3j El [15] . As demon- 
strated here this is a consequence of the ferromagnetic 
coupling between the molecular spin and the magnetic 
substrate. 

The measurements were performed in a STM operating 
below 10 -10 mbar at a temperature of 4.6 K. The mag- 




FIG. 1: CoPc molecules adsorbed on cobalt nanoislands 
grown on Cu(lll). a) Structure model for CoPc (hydro- 
gen: white, carbon: gray, nitrogen: violet, cobalt: blue), b) 
Pseudo three-dimensional STM image (40 x 20 nm 2 , 0.1 V, 
0.5 nA). The spatial oscillations on the Cu(lll) surface are 
due to the scattering of the Shockley surface-state electrons. 
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FIG. 2: Spin-polarized tunneling conductance through the 
nanoislands and CoPc. a) Typical spin-polarized dl/dV over 
two cobalt nanoislands of opposite magnetization (orange: jj, 
brown: jj). Feedback loop opened at 0.6 V and 0.5 nA. 

b) Differential conductance (dl/dV) over the centre of single 
CoPc molecules adsorbed on cobalt nanoislands of opposite 
magnetization (dark blue: jj, light blue: jj). Two sets of 
spectra acquired with distinct tips (noted 1 and 2) are pre- 
sented. The spectra acquired with tip 1 are displaced upward 
by 6 nS for clarity. Feedback loop opened at 0.6 V and 0.5 nA. 

c) Asymmetries arising from opposite magnetizations: CoPc 
(dark blue) and Co nanoislands (orange). The asymmetries 
are an average of all the recorded asymmetries obtained with 
different tips. 



netic substrate was obtained by evaporating cobalt onto 
a pristine Cu(lll) single crystal held at room temper- 
ature, leading, at submonolayer coverage, to the forma- 
tion of triangular-like nanoislands two-atomic layers high 
(Fig.fTJ)). Cobalt islands were chosen since they are a well 
studied system p~6j[17], with, at low temperature, a mag- 
netization perpendicular to the surface. In order to sense 
the spin-dependent contribution through the differential 
conductance, Co-coated metal tips with an out-of-plane 
spin-sensitivity were employed. The dl/dV versus sam- 
ple bias (V) was recorded by superimposing a sinusoidal 
modulation to the junction bias of amplitude 5 mV rms at 
a frequency of 7 kHz, and detecting the first-harmonic of 
the current through a lock-in amplifier. Low-temperature 
constant-current images showed that CoPc adsorbs pref- 
erentially on cobalt, either on top of the nanoislands or 



along their step edges (Fig. [TJ)). Typical island coverage 
was roughly 10 molecules for an area of (10 nm) 2 . On the 
islands, CoPc exhibits a four- lobe pattern 1.4 nm wide 
consistent with the ideal structure of the molecule. 
The spectra of the differential conductance (dl/dV) ver- 
sus sample bias in the center of the islands in the absence 
of CoPc are marked by a dominant SP resonance falling 
at —0.28 V below the Fermi energy (Fig. [2^l), which is be- 
lieved to arise from the hybridization of s — p states with 
the minority d 3z 2_ r 2 band [17 j. Its amplitude is sensitive 
to the spin polarization of the Co islands, so that islands 
with opposite magnetization can be discerned through 
single spectroscopy or through a contrast in spectroscopic 
maps (see Fig. ^jp where V = —0.29 V). In the following, 
islands presenting a strong dl/dV signal at —0.28 V are 
arbitrarily designated as parallel (noted jj), while the 
remaining as antiparallel (jj). All the spectra presented 
were recorded over islands of hep stacking with lateral 
dimensions larger than 12 nm. This ensures that the 
spectroscopic contrast is of SP nature, since structural 
effects as stacking [16 or size-dependent relaxations [17] 
are eliminated. 

Spin-polarized spectra of CoPc were acquired by posi- 
tioning the tip over the Co of the molecule (Fig. [2Jd). A 
broad resonance centered at —0.19 V is detected for all 
CoPc molecules residing on jj islands, whereas a reso- 
nance of weaker amplitude is found for molecules on jj 
islands, indicating that the cobalt atom is spin-polarized. 
Here after, CoPc molecules are noted jj if residing on 
parallel islands, and j j if residing on antiparallel islands. 
The spin polarization of a single CoPc is therefore linked 
to the one of the island, which then indicates that its 
direction may be switched by inverting the island mag- 
netization. This is particularly appealing in view of spin- 
dependent molecular electronics as suggested recently by 
a study on ultrathin films of metalloporphyrin [8] . While 
the resonance is fairly reproducible, the structure at lower 
biases was found to be tip dependent (Fig. [2]o) . To get a 
clearer picture, we have plotted the asymmetries of both 
CoPc (blue line on Fig. [2^) and the cobalt islands (orange 
line) averaged over 13 different tips. The asymmetry is 
defined as (jj — jj)/(jj + jj), the arrows referring to a 
dl/dV acquired on jj and jj CoPc molecules or islands, 
respectively. An oscillatory behavior is observed with a 
sign reversal occurring at —0.26 and —0.64 V for CoPc, 
and at -0.12, -0.42 and -0.80 V for the islands. 

A visual rendering of the dl/dV spectra is presented 
on Fig. [3J where the pixel intensity on each image cor- 
responds to the differential conductance at a given sam- 
ple bias. These images, or dl/dV maps, show isolated 
CoPc molecules on two islands of opposite contrast. A 
marked island contrast (areas without CoPc) is visible in 
Fig. [3^i since the bias is then —0.29 V which corresponds 
to a high magnetic asymmetry (see Fig. [2^). The right 
island is jj, whereas the left one is jj. A similar con- 
trast is observed in the absence of CoPc (Fig. |3b). On 
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FIG. 3: Spin-polarized tunneling conductance maps for CoPc. 
Maps a), c), e) were taken at —0.29, —0.16, and —0.32 V re- 
spectively (image size: 40 x 20 nm 2 ). Each map is normalised 
to span the same colour palette range. Feedback loop opened 
at 0.6 V and 0.5 nA. b) Spin-polarized map of Co nanoislands 
at —0.29 V - as in a) - in the absence of CoPc molecules 
(45 x 25 nm 2 ). d) Profiles of the differential conductance 
along one axis of CoPc extracted from c) (dark blue: |f, light 
blue: ||). f) Map of a f j molecule of e) with model molecular 
structure superimposed (2.6 x 2.6 nm 2 ). 

the contrary, CoPc molecules have similar intensities on 
both islands, as expected by the molecular asymmetry 
which is close to zero at this bias. A contrast between j j 
and t| molecules becomes visible when moving upward 
or downward in energy, as one may expect from the be- 
havior of the asymmetry depicted in Fig. While the 
islands have nearly same intensity at —0.16 V (Fig. [3^), 
the CoPc asymmetry is now maximal. The molecules 
appear as round dots, with a higher intensity for CoPc 
molecules residing on the j j island compared to CoPc on 
the t| island (Fig. [3]i). A contrast is also detected at 
—0.32 V (Fig. [3^), but is now reversed compared to the 
previous map. The benzene ring of the organic ligand 
(Fig. [3]f) appears here as a dim extension of a bright area 
centered at the Co atom of full width at half maximum of 
0.5 nm. We observed that the SP signal, when present, is 
spatially limited to this area (see Fig. SI in the auxiliary 
information) . 

In order to support the experimental findings of a 
twofold spin orientation for CoPc and gain some insight 
on the magnetic coupling with the island, we carried out 
first-principles calculations based on density functional 
theory by means of the PWSCF package [18] . A detailed 
presentation of the method is provided as auxiliary in- 
formation. The magnetic surface is mimicked by using 
a slab model consisting of three Cu layers with two Co 
layers on top and about 2.2 nm vacuum to the next pe- 
riodic repeated layer. The CoPc was placed above the 
Co-layers resulting in a 302 atom model. During the ge- 
ometry optimization only the Cu atoms where held fixed 



at their positions according to bulk values, whereas CoPc 
and the two Co-layers where allowed to fully relax. Dif- 
ferent starting geometries with the cobalt atom of CoPc 
in a top, bridge and hollow position have been checked. 
The lowest energy from the geometry optimization has 
been obtained for the bridge position. The calculated 
distance between the surface atoms and the cobalt atom 
of CoPc is about 0.25 nm, close to the Co-Co distance 
in the surface layers, with only a small distortion of the 
surface occurring (Fig. S2). 

The spin-polarized partial density of states (PDOS) 
gives information on possible hybridizations among 
atoms of CoPc and with the surface atoms (Fig. |4^i). 
The broad DOS between —4 eV and the Fermi level 
(Ep) is dominated by the hybridization of the d n orbital 
(dxz,dyz) with N p^-states (explicitly shown by a green 
solid line on Fig.[4]o). The remaining states are realized 
by d x 2_ y 2- and d^-states (dark red line on Fig. [4)3). The 
prominent minority peak at — 2 eV corresponds to d xy . A 
bonding combination of p-N and d x i_ y i falls at —6 eV for 
majority states (Fig. [4^). A small hybridization between 
cobalt surface atoms and nitrogen/carbon atoms of CoPc 
providing the chemical bonding to the surface is found at 
binding energies between —4 to —5 eV. Due to this in- 
teraction the CoPc is no longer planar (Fig. S2). The 
d z 2 state shows strong hybridization with p-N which dis- 
tributes it over several occupied states (dark blue line on 
Fig. ^p). In contrast to a free CoPc molecule, where the 
spin down lowest unoccupied molecular orbital (LUMO) 
can be identified as d z 2 [T2] . it becomes occupied here. 
This results in a reduced magnetic moment of about 
0.7 jiB for CoPc adsorbed on the cobalt island. In order 
to obtain a more detailed insight into the chemical bond- 
ing one can analyze the wavefunctions or charge densities 
within an energy window, which reveal that there is also 
direct overlap between the surface Co surface atoms and 
the Co atom of CoPc, as shown in Fig. S3. The PDOS 
shows in particular that the resonance at —0.19 V de- 
tected by SP-STS over the cobalt atom of CoPc likely 
originates from a mixture of minority d z i and d n states 
(Fig. [IJd), the contribution of d x i_ y i and d xy - states be- 
ing negligible near Ep. This conclusion should also hold 
in the vacuum region where the tip is, with however a 
possible change in the relative amplitudes of the d z i con- 
tribution with respect to the d n one, since states with 
a small in-plane k -component are the states that show 
the slowest decay into vacuum, i. e. the states to which 
STM is the most sensitive. 

The magnetization density (Fig. [4^) shows clear ferro- 
magnetic coupling between CoPc and the magnetic sur- 
face. While the largest magnetization density is found 
close to the cobalt atoms, as expected, there is also some 
very small magnetization density at some nitrogen and 
carbon atoms (< 0.05 which might explain why a 
SP signal could not be evidenced over these atoms. Some 
of the outer ligand atoms of CoPc show antiparallel align- 
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FIG. 4: Electronic and magnetic structure of CoPc on the 
cobalt nanoislands. a) Spin-resolved density of states (DOS) 
of the Co atom of CoPc (orange line), of a nitrogen atom 
(black line), and of a cobalt atom of the surface with the 
Co-N-CoPc angle close to 89.5° (red line). Inset: Schematic 
adsorption geometry of CoPc. b) Spin-resolved DOS of the Co 
atom of CoPc and contribution from the d- states involved, c) 
Isosurface plot of the magnetization density (red colour corre- 
sponds to spin up and blue spin down). The CoPc is coupled 
ferromagnetically to the surface. There is a spin polarization 
on some carbon and nitrogen atoms. 



merit of the magnetization, which is characteristic for 
spin polarization through these atoms. The bridge po- 
sition is driven by the interaction of the four N atoms 
around the central cobalt atom of CoPc (inset of Fig. [4^i 
and Fig. S2). The geometry optimization reveals that one 
of the N atoms is on top of a surface atom. The remain- 
ing three N atoms are in less favorable positions due to 
the misfit of the molecular geometry and the cobalt sur- 
face below. The N atom on top of a surface atom should 
show the largest wavefuction overlap. The angle between 
the Co surface atom, the N atom on top and the cobalt 
atom from CoPc is 89.5°, which would favor ferromag- 
netic exchange according to the Goodenough-Kanamori 
rules for superexchange [19] . This indirect exchange cou- 
pling mechanism is in accordance with the conclusions 
drawn for metalloporphyrins JH], here, however, the di- 
rect Co-CoPc distance is comparable to the average dis- 
tance within the cobalt surface, so that direct exchange 



may contribute as well but in a weaker extent. 

With the ultimate lateral resolution of SP-STM we 
have shown that the spin state of a single molecule can 
be directly visualized by means of spin-polarized tunnel 
electrons. This is possible because a spin-polarized reso- 
nance built on d molecular orbitals is found close to the 
Fermi level. Additionally the orientation of the molec- 
ular spin is fixed in time owing to a ferromagnetic cou- 
pling with the magnetic substrate. The organic ligands 
play an important role in both aspects, suggesting that 
control over spin transport may be achieved through co- 
ordination chemistry. Our results also open the excit- 
ing perspective of observing changes in single-molecule 
magnetism induced by external stimuli such as light or a 
magnetic field. 
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